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Abstract 

We present an approach for accurate calculation of exchange en- 
ergy and total energy for excited states using time independent density 
i-£h \ functional formalism. This is done by inclusion of gradient correction 

into the excited state exchange energy functionals developed by us. 
We have incorporated Becke and Perdew Wang corrections into our 
functional and have studied various types of excited states having one 
and two gaps in occupation of orbitals. 

1 introduction 

>v 

Exchange energy functional for excited states which we reported in our ear- 
lier works [1] does give accurate excitation energy but absolute values of 
total energy it leads to is not as accurate. Therefore, further correction 
needs to be made in this functional so as to obtain correct total energy of 
the sytem. Inclusion of gradient corrections to improve the accuracy of en- 
ergy functionals has been in practice since the inception of density functional 
theory [2j [3j E]. In absence of exact energy functional for interacting elec- 
tron systems, this is the only option one is left with. Such corrections are 
either derived from the properties of many electrons system or constructed 
impirically or semi-impirically. Generalized gradient approximation(GGA) 
and meta-generalized gradient approximation(meta-GGA) respectively are 
the examples of such corrections |2j [3]. While GGA uses gradient of density 
and its higher power, meta-GGA includes laplacian of density and its higer 
power. Perdew Wang [3] and Becke's j2] correction to exchange energy are the 
example of GGA and Perdew Burke Ernzerhof's [I] correction to exchange 
energy is a meta-GGA kind of correction. 

In Kohn-Sham formulation of density functional theory [7] kinetic energy 
part is well taken care of. And correlation energy which does play important 
role in many instances is only a meager part. Therefore, it is basically the 
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ecxhange energy which needs to be treated accurately at lease so far as we 
are concerned with calculation of total energy of the system. For ground 
states, the local density approximation for exchange energy functional gives 
an everage error of about 10 %. To minimize this error gradient correction for 
exchange has to be used. Among others, the widely used gradient corrections 
for exchange energy are due to Becke [2] and Perdew Wang [3J. Becke's 
correction is a semiempirical one and the PW's correction is obtained by 
gradient expansion of exchange hole. 

2 Gradient correction for excited state ex- 
change energy functional 

Exchange energy functional reported in reference [T] is quite accurate as we 
have seen in the results for excitation energies obtained by using it. To make 
it more accurate we are going to incorporate gradient corrections into it. 
Since doing gradient correction for excited is not as simple as for the ground 
state cases. Therefore we resort to some sort of approximation. We see that 
Becke's and PW's corrections to exchange energy have been constructed to 
satisfy certain properites which are not specfic to the ground states. Such as 
the asymptotic decay of exchange energy density should go as y^. Therefore, 
we hope that these gradient correction will work equally well for excited states 
also. With this intuition we use gradient corrections for exchange energy due 
to Becke ?? and Perdew Wang(PW) ?? to calculate exchange energy and 
total energy for excited states.. 

Exchange energy functional for excited states with Becke's gradient cor- 
rection is given as, 



where, x a = and j3 = 0.0042 and exchange energy functional for excited 





states with PW's gradient correction is given as 





where, 
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Table 1: Total energy and excitation energy of various atoms/ions. 



atoms/ions 

N{2s 2 2p 6 4 S 2s l 2p 4 4 P) 
+ (2s 2 2p 3 4 S -> 2s 1 2p 4 4 P) 
0{2s 2 2p A 3 P -> 2s 1 2p 5 3 P) 
F + {2s 2 2p A 3 P ->■ 2s 1 2p 5 3 P) 
F(2s 2 2p 5 2 P -> 2s 1 2j9 6 2 S) 
iVe+(2s 2 2p 5 2 P -> 2s 1 2j9 6 2 S) 



77 1 * 

~ HF 


77 1 * 

^Becke 




AE HF 


^ E Becke 




53.988 


53.989 


54.026 


0.413 


0.409 


0.423 


73.820 


73.797 


73.835 


0.553 


0.564 


0.582 


74.184 


74.181 


74.227 


0.625 


0.631 


0.652 


98.033 


98.006 


98.051 


0.799 


0.802 


0.0.840 


98.531 


98.561 


98.610 


0.878 


0.869 


0.898 


126.861 


126.745 


126.791 


1.083 


1.074 


1.111 



- —Jk 2 3 - klf In ( ^±^) - —(kl - kf) 2 In ( ^P)(3) 
8tt 3V 3 21 \h-k 2 J 8tt 3V 2 l> \k 2 -k 1 J KJ 

and, s = F = (1 + 0.086s 2 /m + bs A + cs 6 ) m with m = 1/15, b = 14, c = 
0.2 



3 Results 

We have calculated total energies and excitation energies of excited states 
of various atoms using these gradient corrections for exchange enrgies. The 
results obtained in this way are very encouraging, giving total energy for 
excited states which are very close to the corresponging exact exchange only 
values. In the tables [T] to |6] we present total energies and excitation of excited 
states obtained with Becke and PW corrections along with the exact exchange 
only (Hartree-Fock(HF)) values. 

To calculate exchange energy and total energy we do Kohn Sham calcu- 
lation with LSD approximation for exchange potential. From the orbitals 
generated in this way we calculate the k wave vectors from which we get ex- 
chage energy. We use this exchange energy to get total energy for the atoms 
or ions. 



4 Two gap systems 

We have applied the gradient corrected exchange energy functional to excited 
states with two-gaps also. In these cases we use the exchage energy functional 
for two-gap systems has the same form as the one-gap case but has more 
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Table 2: Total energy and excitation enery of various atoms/ions. 



atoms/ions -E* HF - E* Becke -E* PW AE HF AE Becke AE PW 

Li{2s l 2 S 2p x 2 P) 7.365 7.357 7.366 0.068 0.0702 0.075 

Na(3s 1 2 S ->■ 3P 1 2 P) 161.786 161.804 161.890 0.073 0.0763 0.083 

Mg + {3s l2 S ^3p 12 P) 199.214 199.210 199.297 0.158 0.171 0.182 



Table 3: Total energy and excitation enery of various atoms/ions. 



atoms/ions 




77>* 

~ HF 


rp* 

^Becke 


TP* 

-^pw 


AE HF 


AE Becke 


AEpw 


P{3s 2 3p i A S -r 3$^ 4 J 
S{3s 2 3p A 3 P -)> 3,8^ 3 F 


') 


340.417 


340.400 


340.507 


0.302 


0.307 


0.324 


') 


397.079 


397.053 


397.165 


0.426 


0.436 


0.454 


Cl+(3s 2 3p 4 3 P 33^ 


3 P) 


458.523 


458.484 


458.596 


0.526 


0.544 


0.564 


Cl(3s 2 3p 5 2 P 33^ 2 , 


S) 


458.917 


458.900 


459.015 


0.565 


0.567 


0.585 


Ar+(3s 2 3p 5 2 P -> Ss^p 6 


2 S) 


525.598 


525.574 


525.686 


0.677 


0.671 


0.703 



Table 4: Total energy and excitation enery of various atoms/ions. 

atoms/ions -E* HF - E* Becke -E* PW AE HF 

P(2s 2 3p 3 4 S 2s 1 3p 4 4 P) 333.834 333.708 333.766 6.882 

S(2s 2 3p 43 P ^2s 1 3p 53 P) 389.257 389.116 389.174 8.246 

Cl + {2s 2 3p A 3 P -> 2s 1 3p 5 3 P) 449.234 449.079 449.132 9.812 

Cl(2s 2 3p 5 2 P -)> 2s l 3p & 2 S) 449.768 449.602 449.659 9.714 

Ar+(2s 2 3p 5 2 P -> 2s 1 3p 6 2 S) 514.882 514.699 514.751 11.393 



A E Becke 


AEpiy 


7.000 


7.065 


8.373 


8.445 


9.949 


10.028 


9.864 


9.941 


11.566 


11.638 



Table 5: Total energy and excitation enery of various atoms/ions, 
atoms/ions 

B(2s 2 2p l 2 P =► 2s x 2p l 2 D) 
C + (2s 2 2p 1 2 P 2s 1 2p 2 2 P) 
C(2s 2 2p 2 3 P -> 2s 1 2p 3 3 D) 
N+(2s 2 2p 2 3 P 2s 1 2p 3 3 P>) 
5i+(3s 2 3p 1 2 P ->■ Ss^p 2 2 P) 
5i(3s 2 3p 2 3 P -> Ss^p 3 3 P) 



Z7 1 * 
~ HF 


^Becke 




AE HF 


AE Becke 


AEpw 


24.312 


24.304 


24.330 


0.217 


0.211 


0.220 


36.963 


36.945 


36.975 


0.329 


0.328 


0.342 


37.395 


37.376 


37.409 


0.294 


0.303 


0.315 


53.917 


53.434 


53.471 


0.414 


0.439 


0.454 


288.299 


288.289 


288.394 


0.274 


0.282 


0.297 


288.62 


288.606 


288.711 


0.234 


0.246 


0.260 



4 



Table 6: Total energy and excitation enery of various atoms/ions. 



atoms/ions -E* HF - E* Becke -E* PW AE HF AE Becke AE PW 



Be(2s 2 1 S -> 2p 2 l D) 
B{2s 2 2p x 2 P -> 2p 3 2 D) 


14.301 


14.291 


14.300 


0.272 


0.264 


0.288 


24.059 


24.023 


24.037 


0.470 


0.491 


0.513 


C + (2s 2 2p 1 2 P -)• 2p 3 2 D) 
C{2s 2 2p 2 3 P -> 2p 4 3 P) 


36.595 


36.540 


36.552 


0.697 


0.733 


0.765 


36.946 


36.934 


36.948 


0.743 


0.746 


0.776 


iV+(2s 2 2p 2 3 P -> 2p 4 3 P) 
iV(2s 2 2p 3 4 S -> 2p 5 2 P) 


52.865 


52.841 


52.892 


1.023 


1.032 


1.173 


53.222 


53.202 


53.224 


1.179 


1.196 


1.225 


0+(2s 2 2p 3 4 S -> 2p 5 2 P) 


72.829 


72.795 


72.811 


1.544 


1.566 


1.606 


0(2s 2 2p 4 3 P -> 2p 6 4 5) 


73.306 


73.319 


73.339 


1.5032 


1.494 


1.540 


F+(2s 2 2p 4 3 P 2p 6 x 5) 


96.934 


96.935 


96.947 


1.898 


1.883 


1.944 


M#(3s 2 X S -> 3p 2 4 D) 
S(3s 2 3p 4 3 P 3p 6 X S) 
P(3s 2 3p 3 4 S ->■ 3p 5 2 P) 


199.357 


199.362 


199.480 


0.258 


0.257 


0.255 


396.478 


96.455 


396.552 


1.027 


1.034 


1.067 


339.865 


339.834 


339.930 


0.854 


0.874 


0.901 


5i+(3s 2 3j9 1 2 P ->■ 3p 3 2 D) 


287.987 


287.944 


288.034 


0.586 


0.627 


0.657 


Si(3s 2 3p 2 3 P -> 3p 4 3 P) 


288.268 


288.250 


288.342 


0.586 


0.602 


0.629 


C/+(3s 2 3p 2 3 P -> 3p 4 3 P) 


457.795 


457.767 


457.860 


1.254 


1.261 


1.300 



terms. For details the readers may refer [TJ. The results calculated for total 
energy and excitation energies with Becke and PW corrections are given in 
tables |7] along with HF values. 

5 Concluding remarks and future prospects 

We have added gradient correction to exchange energy functional to ob- 
tain better approximation for exchange energy functional. The correction 
term used here have been borrowed from ground state DFT on some physi- 
cal grounds and are found to lead the result in right direction as expected. 
Becke and PW corrections which we have employed here are found to con- 
sistently improve the results for excited states exchange energy, total energy 
and excitation energy. The results obtained from such calculation compare 
very closely with Hartree Fock results. 
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Table 7: Total energy and excitation enery of various atoms/ions having two 
gaps. 



atoms/ions -E* HF - E* Becke -E* PW AE HF AE Becke AE 



PW 



B(2p 1 3p 2 Ml — 3, Ms — 1/2) 23.516 23.441 23.442 1.010 " 1.064 1.108 

N(2p A 3p 1 M L = 2,M s = 3/2) 52.958 52.874 52.884 1.438 1.524 1.565 

F(2p 6 3p 1 M L — 1, Ms — 1/2) 97.034 96.959 96.971 2.371 2.471 2.527 

Ne(2p 6 3p 2 M L = l,M s = 1) 124.594 124.420 124.424 3.948 4.167 4.248 

Na(2s 1 2p 6 3p 2 M L = 1, M s = 1/2) 159.236 159.159 159.203 2.615 2.722 2.770 

Al(2s 1 2p 6 3s 1 3p 3 Ml = 0, Ms = 1/2) 237.16 237.064 237.115 4.702 4.816 4.879 

Ar(2p 6 3s°3p 6 Ap 2 M L = 2, M s = 0) 524.215 524.074 524.156 2.589 2.722 2.776 

K(2p 6 3s 1 3pHp 2 M L = I, M s = 1/2) 597.523 597.436 597.542 1.626 1.850 2.169 
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